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Abstract 

We study the fractional quantum Hall effect in three dimensional systems 
consisting of infinitely many stacked two dimensional electron gases placed in 
transverse magnetic fields. This limit introduces new features into the bulk 
physics such as quasiparticles with non-trivial internal structure, irrational 
braiding phases, and the necessity of a boundary hierarchy construction for 
inter layer correlated states. The bulk states host a family of surface phases 
obtained by hybridizing the edge states in each layer. We analyze the surface 
conduction in these phases by means of sum rule and renormalization group 
arguments and by explicit computations at weak tunneling in the presence 
of disorder. We find that in cases where the interlayer electron tunneling 
is not relevant in the clean limit, the surface phases are chiral semi-metals 
that conduct only in the presence of disorder or at finite temperature. We 
show that this class of problems which are naturally formulated as interacting 
bosonic theories can be fermionized by a general technique that could prove 
useful in the solution of such "one and a half" dimensional problems. 



Typeset using REVTeX 



I. INTRODUCTION 



The quantum Hall effect arises from physics that is special to two dimensions. Neverthe- 
less, it was demonstrated in an early experiment |1| that it survives a small amount of three 
dimensionality; more precisely, that an infinite layer system with interlayer tunneling weak 
compared to the single layer (mobility) gap would continue to exhibit dissipationless trans- 
port and a quantized Hall resistance per layer. The quantized Hall phases in the organics 
also rely on the same effect 0. In the presence of tunneling the chiral edge states which 
exist in each layer hybridize and yield a family of "one and a half" dimensional phases that 
live on the surfaces of the three dimensional systems and exhibit interesting transport in the 
direction transverse to the layers. 

Following the pioneering theoretical work a flurry of work M has focused on the 

integer Hall effect in infinite layers systems with a particular emphasis on the properties 
of the chiral metal that forms at their surface in the quantum Hall phases and some of 
this has found support in experiments ||. Recently, interaction effects and magnetic order 
at v — 1 (per layer) have been discussed as well [j/j. All of this work has antecedents in 
work on multi-component systems such as double layer devices ||, but the limit of infinitely 
many layers sharpens quantitative differences into qualitatively new features, e.g., a different 
universality class for the transitions out of the quantum Hall phases 0. 

In this paper we extend the study of three dimensional quantum Hall phases in two 
directions. First, we offer a systematic analysis of the simplest fractional quantum Hall 
phases in infinite layer systems with a special focus on the phases that exhibit interlayer 
correlations. Here we find several novel bulk properties, most notably a non-trivial structure 
for the quasiparticles. This part of our work builds on the pioneering and early work of Qiu, 
Joynt and MacDonald (QJM) |§, who studied the energetics of multilayer states, and noted 
the irrational partitioning of the quasiparticle charge. The second axis of our work is the 
analysis of the "one and a half" dimensional phases that arise at the surfaces of fractional 
quantum Hall phases. In this regard we report a general analysis of the surface conduction 
by combining renormalization group arguments (which have strong consequences for the 
ground state structure in chiral systems) and a conductivity sum rule. We also carry out 
computations of the surface conductivity of the disordered system in a weak tunneling 
expansion for a large class of states. Together, these show that the surfaces of states with 
marginal or irrelevant electron tunneling are semi-metallic, in that they are perfect insulators 
at all frequencies in the clean limit but begin to conduct at finite temperature and disorder. 



Some of the results derived here were announced previously in a companion paper [nj . 

This last theme, of analyzing the surface phases, is interesting from a more theoretical 
standpoint in that the surface phases are "half" a dimension up from one dimension where 
exact solutions are often possible. We have not made very much progress on this front. What 
we have accomplished is to solve the case of bilayers in some generality > ^ n( ^ special 

solutions for three and four layers, and recast the infinite layer problems in algebraic and 
non-standard fermionized formulations that appear potentially fruitful. We report these here 
in the hope that some readers will find them useful in carrying this program to completion. 
We should also note that the construction of more complex bulk states than those considered 
in this paper could well lead to a richer class of such problems. 

The outline of this paper is as follows. We begin in Section |H] with a discussion of the bulk 
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properties of multilayer fractional states, in particular the novel features of their quasihole 
excitations which include a non-trivial internal charge structure and irrational statistics. 
In Section |TTT| we discuss the edge theory of general multilayer states in the presence of 
nearest-neighbor single-electron tunneling. We specialize to the case where the tunneling 
is marginal in Section [IV], first considering the clean limit and then adding disorder and 
interactions. In Section [V] we extend our analysis to states where tunneling is irrelevant. We 
conclude in Section [VT| by presenting a method for fermionizing the multilayer edge theory 
via the addition of auxiliary degrees of freedom. The appendices contain a discussion of an 
exactly soluble model problem with marginal tunneling (App. [Ap, exact solutions for the 
spectra of two edge theories for the special case of four layers (App. |B|), and some technical 
details including Klein factors (App. and the proof of an assertion made in the main 
text (App. [TJ. 



II. BULK PROPERTIES 



Consider a system which consists of N parallel layers of 2DEGs in a strong perpendicular 
magnetic field. We assume there is a confining potential which restricts the electrons in each 
layer to a region with the topology of a disc. The phase diagram of this system was first 
investigated by Qiu, Joynt, and MacDonald (QJM) JJ. At low electron densities they found 
a variety of solid phases, while at higher densities they found the ground state to be an 
incompressible fluid described by the iV-layer generalization of the Laughlin wavefunction: 

N Ni N Ni Nj 

Mi**}) = n n (** - n n nc*« - e-^** M /4 . a) 

i=l a</3 i<j a=l (3=1 

Here Zi a is the coordinate of electron a in layer i, and Ni is the number of electrons in layer 
i. The exponents are specified by a symmetric, N x N matrix K. The diagonal elements 
of this matrix determine the electron correlations within each layer and the off-diagonal 
elements specify the correlations between layers. For the wavefunction to describe electrons 
the diagonal elements of K must be odd integers. The filling factor in layer i is given by 

JV 

n = (2) 

Note that it is possible to construct more complicated multilayer states by beginning with 
the states in Eq. (HD and carrying out a generalization of the single-layer Haldane-Halperin 



hierarchy construction |13| , [14| . For example, at the first level of the hierarchy the effective 
K matrix is given by 

K {1) =K + AP-\ (3) 

where A is a diagonal N x N matrix with elements ±1, and P is a symmetric N x N matrix 
with even integer elements along the diagonal and integer elements everywhere else. The 
sign of the element Ajj determines whether the excitations in layer j are quasiholes (+1) 
or quasielectrons (—1) and the matrix P specifies the (bosonic) multilayer state into which 
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these excitations condense. Using the matrix Km in Eq. (0) gives the filling factors at the 
first level of the hierarchy. 

We will restrict ourselves to the case where the matrix K is tridiagonal 

K {j = m5ij + n(5i d ^i + S iij+1 ), (4) 

where m and n are nonnegative integers. Since we are interested in the limit of large N, 
we will impose periodic boundary conditions along the direction perpendicular to the layers 
by the identification N + 1 = 1. One of the difficulties encountered when working without 
periodic boundary conditions is discussed briefly below. The standard convention is to refer 
to the state whose K matrix is of the form (|j) as the "nrrai" state. Using Eqns. @) and 
(||), the filling factor per layer in the nmn state is 

v= . (5) 

m + 2n 

From this expression we see that there are multiple states with the same filling factor per 
layer. For example, the states 050 and 131 both have filling factor v = 1/5 per layer. QJM 
found that as the interlayer separation d was decreased there is a phase transition between 
the 050 state in which there are no interlayer correlations and the 131 state in which electrons 
in neighboring layers are correlated. 

Note that Eq. fl5|) applies to the case N —>■ oo or to the case of finite N with periodic 
boundary conditions. For a physically realizable multilayer system, i.e., one with N finite but 
without periodic boundary conditions, one finds that the filling factor Uj is not independent 
of the layer index j. If one solves Eq. (Q) for finite N without periodic boundary conditions, 
using the K matrix given in Eq. (HI) one finds 



1 



m + 2n 
where 



l yt + y N+1 -i 



1 + y 



N+l 



(6) 




(7) 

From this form we find that for large N the filling factor near the middle of the multilayer 
stack is given approximately by Eq. (|5]), but exhibits oscillations about this value as one 
approaches the end layers. Such non- uniformities in the electron density would cost elec- 
trostatic energy, and could be mitigated by the formation of quasihole excitations near the 
outermost layers. Therefore, one would expect that the true ground state for a finite mul- 
tilayer with interlayer correlations would be determined by balancing the creation energy 
of quasiholes against the electrostatic energy of the density oscillations. It is interesting 
to note that to construct a state in a finite stack with uniform filling by carrying out the 
hierarchy construction in only the outermost layers one must go to an infinite level in the 
hierarchy. Henceforth we will avoid these complications by working with periodic boundary 
conditions, arguing that in the limit of large iV they can safely be ignored. 
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A. Quasihole Excitations 



In the last section we learned that for a range of electron densities the ground state of 
the multilayer system is an incompressible state described by the wavefunction flU). It is 
natural to ask about the properties of the excitations of this state. By analogy with the 
single-layer quantum Hall effect we write the wavefunction for a single quasihole at point £ 
in layer j as: 

%*)({*«}) = II - O*0({*a})- (8) 

a=l 

If one interprets l^r^nl 2 = e~^ v ^ z ^' as the Boltzmann factor for a classical 2D generalized 
Coulomb plasma at inverse temperature (3 with an impurity at £, the perfect screening 
conditions give: 

K lk Qf = 5a, (9) 

where Q$ = (K~ l ) kj is the local deviation of the charge in layer k from its ground state 
value, due to the quasihole in layer j . The total charge of the quasihole is 

Q = ^ = UK-^ = ^ = ^, do) 

where we have used Eq. (^). Thus, the relation between the total quasihole charge and the 
filling factor is the same as in the single-layer quantum Hall effect. However, if there are 
interlayer correlations, i.e., n^O, the quasihole charge is spread over many layers. Indeed, 
in the limit iV — > oo the individual charges are irrational: 

\k-j\ 



^ k Vm 2 - 4n* ^ 2n ) ' W 

The same Coulomb plasma analogy, combined with the mean field approximation, allows 
us to find the spatial distribution of the extra charge excited by the hole. The Debye 
screening equation for an average "potential" ®(z) acting on an electron at the point z 
in layer i can be written as 

- (5V 2 vl m = ^8{z - - 2 + 4tt J2 K ik n k , (12) 

k 

where the single-particle average charge density 

m = (27rE i ^)~ 1 e~ m (13) 

is chosen so that at Vi = the unperturbed density would be restored. Clearly, the screening 
in each successive layer is limited by the Debye screening length which is on the order of 
the magnetic length (I = 1 in chosen units). Therefore, as the amount of charge induced 
by the hole is reduced with the separation along the stack, this charge also spreads over a 
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wider and wider area, further reducing the charge density. The shape of the induced charge 
distribution can be found by linearizing the density ([13| ) in Eq. fll^) and solving the resulting 
set of linear partial differential equations by Fourier transformation. In particular, for the 
r.m.s. spread of the distribution created in layer j by the quasihole in layer i we obtain 



r 2 



-2h{\j-i\)/h{\j-i\), where I p (j) = fj dq C -^±, (14) 



and K q = m + 2ncosg is the Fourier transform of Eq. (H). Specifically, for the 131 state we 
find (rf) = 4 (3 + y/5\j\)/5. If we were to color in the perturbed charge density, the hole 
would produce a characteristic "hourglass" shape. 

To compute the statistics of these excitations we follow the method of Arovas, Schrieffer, 
and Wilczek [|15| and consider a state with two quasiholes, one at point £ in layer j and one 
at point r) in layer k: 

^(m-av) = n (%<» -on - v)*o- (is) 

a=l 0=1 

The Berry phase for moving the quasihole at (k, rf) around a closed loop of radius R con- 
taining the quasihole at (j, £) is 



7s = i f ^(*0'.C;M)l^7*Cj',e;M)) 



i j> dr] d 2 z -^—(^(j,^k, v )\ J2 S(2) ( z - z w)\^(j&k,v)}- ( 16 ) 

J£ J Tj Z p =1 



The expectation value appearing in this equation is just the electron density at the point z 
in layer k for the two-quasihole state. We can write this as 

N k 

E * (2) (* - * w )l%e*,,)> = pf(*) + rf'^ W + rf°(*), (17) 

where pjj (2) is the density in the ground state, 8pf : ' r> \z) is the change in density due to 
the quasihole at (k, rf) and 8p^\z) is the change in density due to the quasihole at (j, £). 
If we substitute Eq. ( [FTP into Eq. (|H]), the pj^ term gives the Aharonov-Bohm phase, the 
fiPk^ t erm gives zero by symmetry, and hence 

(stat) _ , f « jj a 7" , 1 _ „ n (j) 



=ij<?z 6p%*>(z) jdr^-—^ -2xQ*>, (18) 

where once again is the charge deviation in layer k due to a quasihole in layer j, which 
we found above to be equal to (K' 1 )^, see Eq. @. The statistical angle for interchanging 
two quasiholes in the same layer is 

A7 = \l { r ] = -*(K-% = - * 2 . (19) 

We find the statistical angle of the quasiholes is an irrational multiple of ir in the infinite 
layer system. 
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III. EDGE THEORY 



In this section we consider the edge theory of the multilayer nmn state in the presence of 
interlayer single-electron tunneling. The low-energy effective Hamiltonian of the edge theory 



in the absence of tunneling is [16 



rL/2 1 

H = dx—Vij :d x Uid x Uj:, (20) 

J-L/2 ATX 

where V is a symmetric, positive definite, N x N matrix which depends on the interactions 
and confining potentials at the edge, and we take the x-axis along the edges of length L. The 
N chiral bosons appearing in this Hamiltonian obey the equal-time commutation relations 

[ui(x), Uj(x')] = inKij sgn(x — x'). (21) 

The electron charge density and the electron creation operator in layer i are given by 

Pi (x) = ^-(K-%d xUj (x), *l(x) oc e~ iu ^\ (22) 
Lit 

Note that because of the factor of K~ l in the expression for the density operator Q2"2"D 
there is a non-trivial relation between the matrix V appearing in the Hamiltonian and the 
matrix determining the density-density couplings, which we will call U. If we rewrite the 
Hamiltonian ( p0|) as 

,L/2 

Ho = J dxUij :pi(x)pj(x):, (23) 



which serves as our definition of U, we find by equating these two forms that 

V = -K- l UK-\ (24) 

7T 

We assume that the system is translationally invariant along the direction perpendicular 
to the layers, which we shall take to be the z-axis. We shall often refer to this as the 
"vertical" direction. Specifically, we assume the matrices and Uij depend only on the 
difference \i — j\. The tridiagonal form of K given above in Eq. (f|) certainly obeys this 
constraint, provided we recall that we are assuming periodic boundary conditions along z. 
From Eq. (24) we see that if both K and U are translationally invariant, so is V. Hence 
to diagonalize K and V we perform a discrete Fourier transformation in the z direction, 
defining: 



N 



U q {X) 



Ee-^Cn), (25) 



^=1 



where 

qE {2im/N\n = 0,l,...,N -1} , (26) 
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and we identify q and q + 2irk for any integer k. 

Using transformation (p5|), the commutation relations (^) become 



[u q (x),u q/ (x')} = m6 q+qlfi K q sgn(x - a/), (27) 

where 

K q = m + 2n cos g, (28) 

are the eigenvalues of the matrix (f|). We will largely restrict our analysis to those states for 
which all the edge modes move in the same direction, i.e., the maximally chiral case. This 
requires that K be positive definite. From Eqns. fl26|) and this implies n < m/2. 
Finally, we rescale the fields, defining 

<f> q (x) = -i=>u q (x), (29) 



which have conventionally normalized commutation relations 

[<j> q {x), <j)q'{x')\ = iiiS q+q ^ Q sgn(x - x'). (30) 



Using the transformations (pq) and ( P9|) to express the Hamiltonian (|20|) in terms of the 
fields g we find 

rL/2 I 

H = dx — v q : d x (j) q d x (j)- q : . (31) 

J-L/2 47T 

The velocities v q are given by v q = V q K q where V q are the eigenvalues of the V matrix. We 
now specialize to the case where the density-density coupling matrix is of the form 

Uij = irvSij + ngiSij-i + S i)j+1 ), (32) 

where v is a velocity determined by the confining potential at the edge and g parameterizes 
the nearest-neighbor density-density interactions between the layers. Using the relation 
between U and V ( |24| ) we find that the eigenmode velocities appearing in the Hamiltonian 



(pT|) are 



v + 2gcosq 

v q = . (33) 

m + 2n cos q 



Next we consider adding interlayer single-electron tunneling to the multilayer edge theory. 
The electron annihilation operator in layer i is given in Eq. (p2[) in terms of the original 
bosonic field Ui{x). Using the transformations ( p5|) and (|29| ) the tunneling operator between 
layers j and j + 1 can be written 

X^j(x)¥ j+1 (x) + h.c. oc \e ia ^"^ + h.c, (34) 

where A is the tunneling amplitude and we have defined the coefficients 
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a n = -=^{\-e^K q . (35) 



As defined in Eq. (B2p the electron operators in different layers do not obey proper anti- 
commutation relations. In Appendix we demonstrate that this can be remedied without 
significantly altering the form of the tunneling operator (|34|). 

The lowest-order perturbative RG flow of the tunneling amplitude A, assumed to be 
uniform along the edge, is controlled by the scaling dimension, 5\, of the tunneling operator 

(0): 

§ = (2 - h) A, (36) 

where the short-distance cutoff increases as I increases. Since the fields <p q obey canoni- 
cally normalized commutation relations, (|30|), we can read off the scaling dimension of the 
tunneling operator from Eqns. (p8|) , (|34"D , and ([35]): 



S *=oJ2 a ii a J-i = m-n, (37) 

Z q 

where there is no sum on j. Using Eqns. (|3Bp and Q3"7| ) we see that tunneling is relevant for 
m < n + 2, irrelevant for m > n + 2, and marginal for m = n + 2. 

If we combine this result with the condition of maximum chirality, n < m/2, we obtain the 
diagram in Fig. |l| The only maximally- chiral state for which the tunneling is relevant is 010, 
i.e., uncorrelated v = 1 layers. There are two maximally-chiral states for which tunneling 
is marginal: 131 and the bosonic state 020. For all other maximally-chiral multilayer states 
interlayer electron tunneling is irrelevant. 

There are two bosonic states, 121 with relevant tunneling and 242 with marginal tun- 
neling, which lie on the boundary of the maximally-chiral region. The corresponding K 
matrices have zero eigenvalues in the limit of an infinite number of layers. The complication 
this zero eigenvalue adds is the possibility of "soft" modes which are not strictly confined 
to the edge Jr7| . The 121 state is the first of a sequence of "pyramid" states: 121, 12321, 
1234321,. . ., where in an obvious extension of the nmn notation the state onmno has next- 
nearest-neighbor correlations specified by the exponent o, in addition to nearest-neighbor 



(n) and intralayer correlations (m), and similarly for the ponmnop states, etc [0. All of the 
pyramid states lie on the boundary of the maximally chiral region (i.e., their K matrices are 
positive semi-definite) and have relevant tunneling. Were it not for the complication from 
the soft modes, their edge theories would be exactly soluble by fermionization as discussed 
in Section |VT]. 
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FIGURES 




FIG. 1. The m-n plane. States to the right of the solid line are maximally chiral. The 
dashed line separates states with irrelevant tunneling (to the right) from states with relevant 
tunneling (to the left). The shaded region contains all maximally-chiral states with non- irrelevant 
tunneling: 010, 020, and 131 (dark circles). The states 121 and 242 which lie on the boundary of 
the maximally-chiral region are also shown (open circles). 



The full Hamiltonian of the multilayer edge theory with tunneling can be written using 
Eqns. ©, ©, and (0): 



n 



L/2 
-L/2 



dx 



1 

47T 



A 



(38) 



where a is a short distance cutoff. In the following sections we will be concerned with the 
z-axis conductivity of this theory, so we will need the form of the current operator along 
z. If pf D \x,t) is the two-dimensional charge density operator in layer j, the continuity 



equation can be written 

: pf D \x,t) = - (Jj |. ; (./'. / S - Jj.j . | (.''. / )j ilJJU.I). 



d_ 

di' 



(39) 



where J* is the current density along the edge of layer j, J* x is the current density flowing 
from layer j to layer I, and d is the layer separation. Introducing the discrete Fourier 
transforms 



p( 2D \k,q,t) = dY^J dxe- iqj e- ikx pY u \x,t), 
j 

with an analogous definition for J~ x (k,q,t), and 

J z (k,q,t) = d^J dxe-^-^e-^J^/x^), 



(40) 



(41) 



where the allowed values of the dimensionless ^-momentum q are given in Eq. (^6|), the 
continuity equation ( |3"5| ) becomes 
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^-/ 2D) (k, g, t) = ~ sm(q/2)J z (k, q, t) - ikj x {k, g, t). 



(42) 



Note pf (x,t) = (l/d)pj(x,t), where pj is given in Eq. fl2"2"|). Therefore, at = we can 
write Eq. @ as 



J z {k = 0,q,t) 



2sin(g/2) 

Evaluating the commutator with the help of Eqns. 



H,p^ D \k = 0,q,t) . 

and (134T), we find 



(43) 



' 5 Vt_3/ 5 



J*(t) = jT(A; = 0,g = 0,t) 



zA<i 



(2na) 



E 



L/2 
-L/2 



e ia jq ,d> q ,(x,t) _ h _ c _ 



(44) 



This expression will be used to calculate the z-axis conductivity of the multilayer system. 



IV. MARGINAL TUNNELING CASES 

We know there are only two maximally-chiral states with marginal tunneling: 131 and 
the bosonic state 020. In this section we begin by showing that in the absence of disorder 
these states exhibit insulating behavior in the vertical direction, i.e., a zz (u) is identically 
zero for all frequencies at T = 0. We will then find that adding disorder to the edge 
theory increases the vertical conductivity, and therefore the surface behaves like a "chiral 
semi-metal". We conclude this section by considering how the transport is modified by 
nearest-neighbor density-density interactions. 



A. Clean Limit 

Our discussion of the system in the absence of disorder will proceed in two steps. First, we 
derive a sum rule which relates the frequency integral of the conductivity to the expectation 
value of the tunneling term in the Hamiltonian. Next we prove that there exists a finite 
range of values for the tunneling amplitude A, and the interaction strength g, for which the 
ground state in the presence of tunneling and interactions is identical to the ground state in 
the absence of tunneling and interactions. We then combine these two results to arrive at 
our conclusions about the zero temperature conductivity. We emphasize that this argument 
is non-perturbative in A and g. This is an important point since we do not have a quadratic 
form of the Hamiltonian for the multilayer state with marginal tunneling. 

The conductivity sum rule is derived by considering the double commutator of the Hamil- 
tonian with the charge density operator. Specifically, we calculate: 

[H,p^ D \k = 0,q,t)},p^ D \k = 0, -q,t)} , (45) 



where the Hamiltonian is given in Eq. (|38|) , and from Eqns. ([22"D , ( p5[) , (29), (flOP, and 
Pj\x,t) = (l/d)pj(x,t) we have 
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1 _ rL/2 1 

(k = 0, q, t) = —VN / dx ^=d x <J) q (x, t) 

27T J -L/2 K 



(46) 



Using the commutation relations for the fields <f> q {x) ( |3"0"| ) it is a straightforward exercise to 
evaluate the commutator (l|5|) with the help of Eqns. fl38|) and (|46|). One finds: 



[H, p {2D) (k = 0, q, t)],p^(k = 0, -q, t) = 4 sin 2 (g/2) Hy 



(47) 



where 7i\ is the tunneling part of the Hamiltonian. 

We next relate the expectation value of this double commutator to the integrated con- 
ductivity. The continuity equation (|39D introduced in the previous section can be written 



9 r dtii r p -iu{t-t') 

(k = 0, q, t) = - sin(g/2) J ^ J dt' e ——J*(k = 0, q, t) 



(48) 



Suppressing the argument k = 0, using [H, p] = —idp/dt, and taking the expectation value 
of Eq. © we find 



[H,pf D) (g,t)],pP(-g,t)l) = -sin 2 (g/2) / =1 / dt e^{[J*{q, 0), J z {-q,t)])- (49) 



d? 



du 



2tt 



If we take the expectation value of Eq. (|4"T|), and use it to eliminate the l.h.s. of Eq. 
find in the limit q — ► 0: 

CHx) = ^ J ^ J dte^{[F(0),F(t)]), 



we 



(50) 



where we have used I z {t) = J z {q = 0,t). The Kubo formula relates the commutator 
appearing in this expression to the z-axis conductivity: 



a zz (cu) 



1 



2u NLd 

Thus we arrive at the final form of the sum rule 



dte ZUJt ([I z (t),r(0)]). 



(51) 



du^ea zz {uj) 



7TC 



NL 



(52) 



This is an exact relation valid at any temperature. It applies to all multilayer states, even 
in the presence of disorder, regardless of the relevancy of the tunneling term. 

Next we turn our attention to proving the stability of the ground state in the presence 
of interlayer tunneling and interactions. We will describe in detail the case of the bosonic 
020 state. A proof along the same lines can be constructed for the 131 state. We begin by 
writing the Hamiltonian of the 020 edge theory in terms of the original radius R = 1 Bose 



fields Uj(x): 



n 



020 



i dx 

J-L/2 V 



167T 



{d x Uj) 2 : +-^-d x Ujd x u j+1 + 



(2na) 



(53) 
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Instead of performing the discrete Fourier transformation along z fl25|) we simply rescale the 
fields, defining Uj(x) = y/2tpj(x), where (fj(x) are radius R = l/y/2 chiral bosons. We can 
then define a set of N su(2)i Kac-Moody (KM) currents: 



J/(s) 



2tt\/2 



d x cpj(x) 



(54) 



in terms of which the Hamiltonian (p3 ) can be written 

i-L/2 



n 



020 



dx 



L/2 



y :[J,(x)] 2 : +ngJ*(x)J* +1 (x) 



+ \{J-{x)Jf +1 {x) + Jf +1 {x)Jf{x)) 
where we have employed the identity 



[L/2 2 rL/2 1 

/ dx :[J z {x)\ := dx — :[J(x 

J-L/2 ' J-L/2 3 



L/2 
-L/2 



(55) 



(56) 



Consider the theory without interactions, g = 0, and without tunneling A = 0. The zero- 
energy ground state of this Hamiltonian, |0), is formed by taking the tensor product of the 
highest-weight state of the spin-0 irreducible representation of each su(2)i algebra. Thus, 
the ground state satisfies 



J? n |0)=0, for Vj; a = x,y,z; n > 0, 
where the Fourier components of the currents are defined by 

pL/2 



J a hn = J L2 dx J*(x)e~ 2mnx/L , neZ, 



and obey the algebra 



n 



[ J i,n> J j,m\ ~ °n+m,0 T iOjjt J j^ n+rn - 



(57) 



(58) 



(59) 



Using Eqns. ( (5"5| ) and ( |58| ) we can write the Hamiltonian as 

A 



H 



020 



TCV 



. ja ja . ,^9_ jz jz | 2_( J- J+ _|_ f~ J+ ) 

6L ' i> n 3~ n ' J' n j^\ u 3,n' J j+\,-n ' u j+l,n u j,-n) 



(60) 



For later reference we also record the expression for the current operator in terms of the KM 
generators 



iXd 



° 3,n° j+l,-n u j+l,n u j,~n 



(61) 



From Eqns. (|57|) and ( |60"D we see that not only is |0) the zero-energy ground state of the 
Hamiltonian with g = A = 0, it is a zero-energy eigenstate of the Hamiltonian for all g and 
A: 
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ft 2o|0> = 0. 



(62) 



Of course, this does not mean that |0) is the ground state of the Hamiltonian with 
non-zero g and A. To establish this we now demonstrate that for a range of g and A the 
Hamiltonian is positive semi-definite. For A > 0, g > 0, we rewrite Eq. (|55|) as 



: J J (x) J+ (x) : + : Jf (x) ./ ; (x) : 



J- 



TTV 

1 T 



: -A 



(63) 



The two terms in the first line are clearly non-negatively defined, while the terms in the 
second line can be shown to be non-negatively defined for v > 2g + 8\/n using Eq. ( |56f) and 
the basic identity 



(64) 



An analogous construction can be done for A < or g < 0, and we conclude that the 
Hamiltonian (1531) is positive semi-definite provided 



2\g\ + 8|A|/tt < v. 



(65) 



This is also the sufficient condition for the stability of the original ground state, as follows 
from our previous result that the Hamiltonian annihilates the original ground state |0). 

Note that the normal ordering does not present a complication because it amounts to 
subtracting off the unit operator times the expectation value of the Hamiltonian in the state 
|0), which is zero. Therefore, since the Hamiltonian is positive semi-definite and we know the 
state |0) is a zero-energy eigenstate, it follows that |0) is the ground state of the Hamiltonian 
provided Eq. d65[) is satisfied. 

We can now combine the sum rule with our knowledge of the exact ground state to say 
something about the conductivity. At zero temperature the expectation value on the r.h.s. 
of the sum rule (|52|) is a ground state expectation value. Since we know the ground state 
|0) is unchanged by tunneling and interactions and (0|7Yi |0) = we can conclude 



0. 



(66) 



The real part of <j z 
Eq. fl66"D we have 



(u) is dissipative and hence it cannot be negative. Therefore, from 



<T ZZ (CU) 



0. 



(67) 



We reiterate that this is an exact statement for the clean 020 multilayer at T = provided 
g and A satisfy Eq. flBSp. As we mentioned above, a similar result can be shown to hold for 
the 131 state, the other maximally-chiral state with marginal tunneling. The condition for 
the stability of the ground state of the 131 edge theory is also given by Eq. (p5|). In the next 
section we will explore what happens when we add disorder to these systems. 
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B. Disordered Case 



In the previous section we found that the clean multilayer with marginal tunneling ex- 
hibits insulating behavior. In this section we study the marginal tunneling case in the 
presence of disorder. We perform a perturbative calculation of the vertical conductivity via 
the Kubo formula, finding that disorder increases the conductivity. 

In Section 111 we discussed the general edge theory of clean multilayer systems. Our first 
task is to add disorder to this model. We will then use the Kubo formula to find the z-axis 
conductivity. 

The disorder we consider is a random scalar potential Vj(x) in each layer j, which couples 
to the Bose fields via the charge density 



L/2 
-L/2 



rL/2 1 

dxVj(x)pj(x) — J dx —Vj(x)Kj k d x Uk(x) 



(68) 



We assume Vj(x) is a Gaussian random variable, uncorrelated between different layers, but 
for now we will not specify how it is correlated within the same layer. 

Using the Fourier transformation fl2~5[ ) and rescaling (^) of the Bose fields we can write 
Eq. pi) as 



L/2 1 1 

dx- =V q d x (f)„q(x), 



(69) 



where 



V q (x) 



(70) 



Recalling the form of the clean multilayer edge theory from Section [TTI|, our full Hamiltonian 



is 



n 



L/2 
-L/2 



dx 



— v a : d x d> a d x d)^ a : + — - 

4tt q rq r q Y ( 2?ra ) 



(e ia ^"^+ h.c. 



+ ^—r=V q d x 4>- q (x) 



(71) 



We move the disorder term into the phase of the tunneling amplitude by shifting the bosons: 

1 



+ 



VqjK q 



dyV q (y). 



(72) 



If we define the phases 



M3i 



dyV q (y), 



(73) 



then with Eq. ([72] ) the Hamiltonian (|7T| ) is 
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n 



L/2 
-L/2 



dx 



—v q :d x (f) q d x (f)_ q : 

47T 



A 



(2na) 



(e^ {x 



)Ja jq cf> q (x) 



+ h.c 



(74) 



and the current operator (fMj) becomes 

iXd ^ r L / 2 



(2vra) : 



E 



dx 



-L/2 



(75) 



We are now in a position to calculate the z-axis conductivity. We will first evaluate 
the Matsubara two-point function of the current operator (|7"5|). Since Eq. (|T4]) is a theory 
of iV interacting chiral bosons, the best we can do is a perturbative calculation in the 
tunneling amplitude A. Although we will work perturbatively in A, we will be treating 
the disorder exactly. We thus expect the C(A 2 ) result to be reliable, since the disorder 
essentially randomizes the phase of the electrons between tunneling events. More formally, 
if the tunneling amplitude is a coordinate-dependent, delta-correlated random variable, then 
the RG flow (|36|) is modified to 



dA 
~d£ 



(3 - 25 X )A, 



(76) 



where A is the variance of the tunneling. Thus for 5\ = 2 we would conclude that the tunnel- 
ing is irrelevant. The combination Xe 1 " 13 ^ that appears in the Hamiltonian (|74]) is in general 
not delta-correlated, as we will discuss below, but since the case of x-independent tunneling 
is marginal we expect any departure from uniformity makes the tunneling irrelevant and 
therefore justifies a perturbative expansion in A. 

Using the expression for the current operator ( |75|) we evaluate the Matsubara two-point 
function 



C(r) = -(T t /*(t)/*(0)> 

d " X2 Jdxdx'J2( T r(e nAx) :e 



L 4 



ia jq 4> q (T,x) . 



x e 



(< 



-ilk(x') . e -ia kq <j> q {0,x') . 



(77) 



where we have normal-ordered the vertex operators. Note that because the current operator 
(f?"5l) has a factor of A out front, there is an explicit factor of A 2 in Eq. (|77|), and therefore 
to find this function to order A 2 we can evaluate the expectation value using the A = 
Hamiltonian. One finds that only the cross terms in the product of the current operators 
give non-vanishing contributions, and only when j = k. The terms have the form 



n 



\i y 2(3iv q ) sinh (j^-(x — x' + iv q r + ia sgn t] 



»2 • 



(78) 



where a 2 



OLj q OLj- q , 



with no sum on j. This is a very complicated function; a great simpli- 



fication occurs if the eigenmode velocities, v q , are identical for all q. From the expression for 
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v q (|33|) we see that this occurs when g = nv/m. For the two multilayer states with marginal 
tunneling, 020 and 131, this condition gives g = and g = v/3, respectively. We will refer to 
the case where v q is independent of q as the "solvable point." For 020 it corresponds to no 
density-density couplings between neighboring layers while for 131 it occurs at a non-zero 
value of the density-density interaction strength. Note that for more general forms of the 
matrices K and U, i.e., not tridiagonal, the solvable point corresponds to U oc K. 

Working at the solvable point, writing v q = vq, and using our knowledge of the scaling 
dimension of the tunneling operator (|37|), 



E 



O" 



E a jq a j-q 



25, 



(79) 



we find upon disorder averaging Eq. (|7 
2d 2 X 2 



CM 



{2(3v y 



^2 J dx dx' 



Wj{x,x') 



sinh — x' + wqt + ia sgn t] 



•l - 



(80) 



where we have defined the function 



Wj(x,x') 



e i[-yj(x)—yj(x')]_ 

If we take the correlation function of the scalar potential to be 



Vj(x)V k (x') = 5 jk Z(x - x'), 



(81) 



(82) 



then using the definition of the phase 7 3 -(x) 
written 



we find that the function Wj(x,x') can be 



Wj(x, x') = exp 



Vq Jx 



-a / dyi 



dy 2 Zfa - y 2 ) 



(83) 



Note that Wj(x, x') = W(x — x') is independent of j and depends only on the magnitude of 
the relative coordinate \x — x'\. Therefore the sum over j in Eq. (|80|) just gives a factor of 
the number of layers, N, and we can replace the integration over x and x' by an integration 
over the average coordinate, which gives a factor of L, and an integration over the relative 
coordinate y = x — x' . If we define the Fourier transform of the function W(x): 



W{k) 



dxe tkx W(x 



then we can write Eq. (|80D as 
C(t) = 



^W(k) j dy 



(2(3v ) 



Jky 



sinh [ik( y + iv ° T + ia sgnr ' 



(84) 



(85) 



The factor of L in this expression will cancel against the factor of 1/L in the Kubo formula 
(pl|), and in anticipation of this we have extended the range of the y integral to include the 
entire real axis since we are interested in the result in the limit where the system size L — > oo. 
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The integration over y can now be performed by the method of residues. The integrand has 
fourth-order poles on the imaginary y-axis at the points y n = i[v (/3n — r) — asgnr] for 
integer n, and is exponentially small in the upper (lower) half plane for k > (k < 0). One 
finds 



dy 



Jky 



TV 

3 



' (3vQk 



+ 



7T 



sinh [ik( y + iv ° T + ia sgnr 

We now Fourier transform Eq. (j35|) with respect to the imaginary time r, 



,f3vok 



(86) 



C(iO r , 



/? dre !U » T C(r 
d 2 X 2 NL 



48vr 4 ^ 



dpW(p/vo) 



p 3 + 4n 2 p/p 2 
iu n +p 



(87) 



where we have changed the integration variable to p = t^A;. Analytically continuing (iu> n 
uj + «0 + ) we find 



C{u) 



dte^([P(t),P(0)]) 



d 2 X 2 NL 
' 48vr 4 ^ 



dp W(p/uo) 



+ Att 2 p/(3 2 
uj + p + i0 + 



Using this result in the Kubo formula ( |5TD with the help of the identity 



-—— = V- T inS(x), 
x ± i0 + x 



(88) 



(89) 



and the observation that W(k) is real since W(x) is an even function of x, we arrive finally 
at 



, X 2 d 



3?e <7 zz (lo 



487r 3 w 4 



(uo 2 + 4tt 2 T 2 )W(uj/vo). 



(90) 



This is our central result. It is the 0(X 2 ) term in the vertical conductivity at the solvable 
point for both the 020 and 131 states at a finite temperature T, including potential disorder 
which enters through the function W(k). The imaginary part of the conductivity can of 
course be found from Eq. ( p0|) by the usual dispersion relations [p~9[] . 



Let us first consider the limit of a clean system, Vj(x) = 0. From Eqns. 



and 



we see that in this case W(x) = 1 and therefore W(k) = 2n5(k). Hence from Eq. (|90| ) we 
find 



Qvq 



T 2 S{uj) 



(91) 



This is consistent with our result in Section |1V A| that the z-axis conductivity in the clean 
system vanishes at T = 0. An interesting aspect of this result is that it is proportional to 
8(u), i.e., the conductivity vanishes at all u> ^ 0, and the DC conductivity is infinite. One 
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question to ask is whether these features are a consequence of the fact that Eq. ( pT|) is only 
the first non-vanishing term in a perturbative expansion in A. 

Although we cannot definitively answer the question about the finiteness of the DC 
conductivity beyond 0(X 2 ), we can prove that at higher orders in A the conductivity cannot 
vanish at all u ^ 0. First we note that in the clean system the current operator commutes 
with the Hamiltonian in the absence of tunneling: [I z ,H,q] = 0. For example, for the 020 
state this can be seen from Eqns. ( |S0"D and (|STD, using the commutation relations for the 
Kac- Moody currents given in Eq. (|59|). Since the current operator commutes with 7i , the 
expectation value of the current-current commutator evaluated using the A = Hamiltonian 
vanishes, and therefore by the Kubo formula ( pH ) the vertical conductivity must vanish for 
all non-zero frequencies. This explains why the conductivity is zero except at u = in 
Eq. (|91]). By the same reasoning we see that if the current operator commutes with the full 
Hamiltonian, [I z ,TC] = 0, then a zz {uS) would be zero for all uj ^ to all orders in A. The 
converse of this is also true, i.e., if a zz [iS) is zero for all u ^ then [I z ,1i] = 0. This is 
proven in Appendix [D[ For the multilayer edge theory with marginal tunneling: [I z , 7i] ^ 0. 
For the 020 case this can be established using Eqns. ( |60| ) and (pip . Therefore, since the 
current operator does not commute with the full Hamiltonian, we can conclude that a zz (u>) 
cannot be zero for all uj ^ for the clean 020 or 131 multilayers. Hence, we would expect 
that the 5(u) factor in Eq. (|91|) would be broadened by the higher-order corrections in A. 
Recalling that A is dimensionless and the conductivity should not depend on the sign of A, 
one plausible scenario is: 

Se a~ {a) = g Tf (w/r) i „ (92 ) 

This conjectured form has several interesting features. First, we find that the DC con- 
ductivity goes to zero with temperature as T. We expect this to be a generic feature of the 
exact result in A for the clean system, provided the DC conductivity is finite. This differs 
from the disordered system ( PO ) where we find the DC conductivity vanishes as T 2 . In 
addition, note that the uj — > limit of Eq. (|9~^) is independent of the tunneling amplitude 
A. 

We now return to the result for a zz (uj) in the presence of disorder fl9~0"|). If the disorder 
is delta-correlated, the function Z(x) defined in Eq. fl82|) is 



Z(x) = A5(x), (93) 

where A is the variance of the disorder. For this case W(x) = exp(— A\x\/vq), which has a 
Fourier transform W(k) = (2A/vg)/(fc 2 + (A/vq) 2 ), and hence Q9"CD yields 

^ , , \ 2 dA uj 2 + 4ti 2 T 2 

a zz (lu) = - = , -= . . . . 94 

Note that at zero temperature the conductivity vanishes in the DC limit. This conclusion 
is independent of the details of the disorder and holds as long as W(0) is finite. We also see 
that the conductivity approaches a constant as uj — > oo. This second feature is a consequence 
of taking the disorder to be delta-correlated. We see that a zz (oj) goes to a constant at large 
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tjj because of the slow fall off of W(k), which is in turn caused by the non-analyticity of 
W(x) at x — 0. We shall now investigate how Eq. ( p4[ ) is modified if the disorder has a finite 
correlation length. 

As a first step, we determine the asymptotic form of the disorder-averaged phase, W(x), 
at small and large x when the scalar potential disorder has a finite correlation length. We 
return to the expression for the correlation function of the random potential (B^) and write 



Z(x) = A 2 A(x/£ ), 



(95) 



where A is a parameter with the dimensions of energy that sets the strength of the disorder, 
Iq is the correlation length of the disorder, and A{z) is a dimensionless function which we 
take to have the properties: A(0) = 1, A(—z) = A(z), and 



rod 

/ dzA(z) 
Jo 



1. 



(96) 



Now consider the disorder-averaged phase, W(x), given in Eq. (|83|). For \x\ -C £ , i.e., 
for distances small compared with the correlation length of the random potential, we can 
replace Z(yi — y 2 ) by its value at y x — y 2 = and obtain: 



W(£ z) 



-(Ae z/v ) 2 



l-^V. for 



\z\ < 1. 



I "0 / 



(97) 



In the opposite limit |x| ^> £q we change integration variables in Eq. ( |33D to y c = (yi + y 2 )/2 
and y r = yi - y 2 : 



W(x) = exp 



Iff f\ x \ 

— 2 / d Vr (\ x \ + yr)Z{y r ) + / dy r 
Vq \J-\x\ Jo 



(\x\ - y r )Z(y r 



(98) 



Since \x\/Iq ^> 1 we can extend the integration over y r to infinity with small error. Using 
the symmetry of A and defining 



we then find 



W(£ z) w exp 



rod 

C = / dz zA(z) 
Jo 



for \z\ 1. 



(99) 



(100) 



We see that at short distances W(x) is parabolic ( p7|) while at long distances it decays 
exponentially ( |100| ). The presence of a finite correlation length for the random potential 
does not change the long distance behavior of W(x), but it does "round out" the non- 
analyticity at x = present in the case of delta-correlated disorder. Generally, this is 
sufficient to make the function rapidly decaying as uj — > oo. 

Having found the generic behavior of the function W(x), we now choose a specific form: 



W(x) 



1 

cosh(£x) 



(101) 
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characterized by the parameter £ which has the dimensions of an energy. This choice of 
W(x) has several appealing features. First, it has the correct asymptotic forms at large and 
small values of x. Indeed, if we match the asymptotics of Eq. ( |101| ) to Eqns. (P?|) and ( |100| ) 
we find 



A 

^0 



1 

2|' 



(102) 



We see that for the choice of W(x) in Eq. ( |101| ), the parameter £ is proportional to the 
strength of the disorder and inversely proportional to the correlation length. Hence this 
form allows us to explore the effect of a finite correlation length in the region A£ ~ 1. 

A second advantage of Eq. (|101|) is that its Fourier transform can be evaluated in closed 
form. A straightforward computation gives 



W(k) 



-ikx 



dx 



7r/£ 



cosh(£x) cosh 



(103) 



Using this result in the general form above for the conductivity ( j90|) we find 

\ 2 d uj 2 + 4tt 2 T 2 



5£e a z 



487r 2 ^£cosh (ttlu/2£v )' 



(104) 



Comparing Eq. ( |104| ) with the expression for the conductivity in the case of delta-correlated 
disorder (|94]) we see that the finite correlation length gives a conductivity which exponentially 
decays to zero as uj —>■ oo rather than approaching a nonzero value. This result for a zz ( |104| ) 
is shown in Fig. [|at various temperatures. For large temperatures the conductivity is peaked 
around u = 0, while at smaller temperatures the maximum occurs at a finite frequency. 



2 . 5 



1.5 




0.5 



1 2 3 4 5 

FIG. 2. The real part of the vertical conductivity for the 020 or 131 multilayer with disorder 
(|104| ). The horizontal axis is frequency measured in units of the parameter vq^. The temperatures 
of the curves, starting with the uppermost, are: T/vqS, = 0.25,0.2,0.15,0.1,0.01. 
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C. Effect of Interactions 



In the preceding section we calculated the order A 2 term in the z-axis conductivity of 
the multilayer states with marginal tunneling, 020 and 131, at the solvable point where 
the eigenmode velocities v q are independent of q. In this section we discuss how the zero- 
temperature result is modified away from the solvable point. 

Since we will only consider T = in this section we shall work directly in real time t. The 
real-time zero-temperature analog of the expression for the Matsubara two-point function 
of the current operator (jSPf) is 



C(t) 



2id 2 \ 2 NL 
(2vr) 4 



dyW(y)H 



1 



(y ~ v qt + sgnt) c 



(105) 



where from the definitions ( [73]) and (|81|) we see that the expression for W(x) in Eq. (|83|) is 
changed to 



W(x) 



exp 



2 ^ sin 
N \ ~ v 2 



f X d yi ( X dy 2 Z{ yi -y 2 ) 
Jo Jo 



(106) 



The corresponding retarded correlation function is 



, 2id 2 X 2 NL 



c R (t) 



(27T) 



0(t) dyW(y) 



n 



n 



(y - v q t + ia) a i q (y - v q t - ia) a i 



(107) 



We would now like to perform the y-integration. From Eq. 
appearing in Eq. (|107Q are 



we see that the exponents 



N 



(3 — cos q — 2 cos 2 q) , 



(108) 



and therefore, as a function of complex y, the first integrand in Eq. (|107| ) has N branch point 
singularities at the points v q t — ia and similarly for the second integrand at v q t + ia. This is 
a complicated singularity structure, but observe that in the first term these branch points 
are all below the real axis, while in the second term they are all above the real axis. Hence, 
since Y^ q a q is an integer, we can certainly choose branch cuts that lie entirely on one side of 
the real axis for each term. When we considered the solvable point in the previous section 
we were able to evaluate the conductivity for any disorder W(y). We shall not attempt the 



same feat away from the solvable point. The evaluation of the y- integral in Eq. ( |107|) is 
greatly simplified if W(y) is a meromorphic function of y. One such function was discussed 
in the previous section (|101|) , and for the remainder of this section we specialize to this form. 

Using W{y) = 1/ cosh(£?/) in Eq. ( |107|) , we can close the first term in the upper-half 
plane and the second term in the lower-half plane, avoiding all the complicated branch point 
singularities, and picking up the residues of the poles of W(y). We find 



CK(t) = - i 6(t)(iP(t),P(0)]) 

^NL em3 - ( 



(2 



7T 



•irn 



~y q &- z?r(r + l/2)] a « 



(109) 
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Therefore, from the Kubo formula 



we arrive at 



\2 j £-3 oo 



(2tt) 3 oj 



[v q ^t-m{r + l/2)}< 



;no) 



This is certainly far from a closed form result, but it is a convenient form for extracting the 
high and low frequency asymptotics of the conductivity. At large (positive) ui the asymptotic 
form is determined by the closest singularity to the real axis in the upper half of the complex 
t plane. From Eq. (|110|) this clearly occurs a finite distance away from the real-t axis and 
we find 



'Stea zz (u) — 5 1 exp 



TYUJ 



111^ 



where v max is the maximum of v q over q. For the 020 state we see from the expression for v q 
(|33|) that at the solvable point (g = 0) Vq = v/2 while away from the solvable point (g > 0) 
Vx = v/2 + g. Comparing Eq. ( p. 11] ) with the asymptotic form of Eq. ( |104| ): 



5£e a z 



exp 



7TUJ 



2£v 



;ii2) 



we find that one effect of interactions at T = is to soften the exponential decay of the 
conductivity at large frequencies. 

Returning to Eq. ( [110] ), we now consider the behavior at small frequencies. If we assume 
to > we can close the t integral in the upper-half plane because the integrand is exponen- 
tially small there. Hence all the terms with r < give no contribution. For each r > there 
remains iV branch point singularities in the upper-half plane enclosed by the contour. Since 
we are interested in small u we stretch the integration contour into a large circle so that 
everywhere along the contour \t\ is very large. We can then expand the integrand in powers 
of 1/t and integrate term- by-term: 



dte iuJt H 



1 



v q it - i7r(r + 1/2)] Q ^ 



9 Vq 1 



i + 



m{r + 1/2) ^ a 



q V 1 



t 2 




1 - 



7T UJ 



+ l/2)E^ + 0( 



(113) 



If we retain only the lowest term in uo in this expansion, we note that upon substituting this 
into Eq. ( |110| ) we would encounter the divergent sum 



r=0 



By considering the solvable point, where v q = v independent of q, it is clear that the proper 
regularization of this sum is 
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y (-1Y = limV (-e~ s Y = lim — x = (115) 



r=0 r=0 

and thus from Eqns. (|110| ), (|113|) , and ( |115|) we have 




ne^(u) ^ — — TT (H6) 



Comparing this with the small frequency asymptotics of Eq. ( |104j ) at T = 

s w _> n X 2 d 1 lu 2 



Me^u) ^ i^-^r, (117) 



we see that the interactions do not modify the result that the conductivity goes to zero in 
the DC limit at zero temperature. While the power of u is the same in Eqns. ( |116| ) and 
( |117|) , the interactions do modify the prefactor. 

To summarize, we have considered how interactions away from the solvable point modify 
the vertical conductivity at T = 0. It should be noted that in arriving at Eq. (|110|) the 
interactions were treated exactly. We find that at large frequencies the interactions soften 
the decay of the conductivity, but it remains exponential, while at small frequencies the 
interactions do not change the frequency exponent of the conductivity. We remark that the 
difficulty encountered when attempting to apply the method presented in this section to 
investigate the effect of interactions at a finite temperature is that in the finite-T version of 
Eq. (|107|) there are complicated branch point singularities on both sides of the real-?/ axis, 
see Eq. ([78]). Nevertheless, we believe that even at a finite temperature the asymptotics of 
<j zz (uj) are not modified by interactions. 



V. IRRELEVANT TUNNELING 

In the previous section we have been concerned exclusively with the case of marginal 
tunneling, namely the 020 and 131 multilayer states. The calculation presented in Sec- 
tion [IV B| can readily be extended to the more general case of the nmn state. Essentially 
the only change in the computation is that the scaling dimension of the tunneling operator, 
5\ = m — n, differs from its marginal value of 2. Therefore, the poles in the complex y-plane 
for the integrand in Eq. (|86| ) are now of order 2(m — n). One finds that for the solvable 
point g = nv/m: 

For the case of relevant tunneling, 8\ — m — n — 1, the product over k is absent in 
Eq. ( |118|) , and if we use W(uj/vq) = (2A)/(u 2 + (A/v ) 2 ) appropriate for delta-correlated 
disorder we reproduce the behavior of the chiral metal. 

In the remaining cases, 5\ > 2, we see at T = the conductivity vanishes as u — > as 
cu 2 ^ Sx ~ 1 \ and at uj — it vanishes with temperature as T 2 ( 5x ~ 1 \ In the absence of interlayer 



24 



correlations, i.e., for the OmO state, this means the DC conductivity obeys a zz ~ T 2m ~ 2 . 
This disagrees with the result of Balents and Fisher, who find by a scaling argument that 
u zz ~ T 2m ~ 3 M. While it is unclear whether Balents and Fisher have in mind the clean 
system or the disordered system, we believe the reconciliation of this discrepancy is that, as 
we found for the marginal case, the temperature scaling of the DC conductivity is different 
in the clean and disordered systems. If we consider the clean limit of Eq. (|118|) by writing 
W{uj/vq) = (vo/T)5(uj/T), and then conjecture that the interactions broaden the delta 
function, we find: 

Ke^M ex ^ ( „ /r) 2 + 0( g2) II + (2^T) 2 ], (119) 

where g is the dimensionless interaction strength. The broadening of the delta function by 
interactions is suggested by the fact that the nearest-neighbor density- density interaction 
Hamiltonian does not commute with the current operator I z . For the case of the OmO state 
the DC conductivity of Eq. ( |119]) goes to zero as T 2m_3 , in agreement with Balents and 
Fisher. 

Finally, returning to the disordered case ( |118| ), we see that for the 050 state the DC 
conductivity scales as T 8 , while for 131 it scales as T 2 Q9"0|). Recall from Eq. (|5|) that both of 
these states occur at the same electron density, v = 1/5 per layer. In their numerical work, 
QJM found a phase transition between these states as the layer separation d was varied. 
The significant difference in the temperature scaling of the DC vertical conductivity in these 
states suggests a way to experimentally detect this phase transition. 



VI. FERMIONIZATION 

In certain special cases the multilayer edge theory can be solved exactly. The case of 
N = 2 layers is discussed extensively elsewhere [jTT| , [12]| , and examples of exact solutions for 
N = 4 layers are given in Appendix [FJ. A potentially useful first step to solving the edge 
theory for arbitrary N is to find a fermionic representation. In this section we discuss the 
fermionization of the edge theory for an arbitrary number of layers N. In some cases the 



procedure involves the introduction of auxiliary degrees of freedom as discussed in Ref. |TT 

We begin with the observation that the scaling dimension ([57]) of the tunneling operator is 
always an integer. This underlies the fact that the chiral Hamiltonian (^) can be fermionized 
exactly, with the tunneling term expressed as a product of Fermi operators whose number 
equals twice the scaling dimension. 

The usual mapping between chiral radius-one bosons and fermions requires independent 
bosonic fields. Even though we are interested in correlated states with non-trivial commu- 
tators (pl|), these can be readily diagonalized by a linear transformation. For example, for 
the 131 state we can write 



it 



j 



fj-1/2 + <Pj + <Pj+i/2, j = l,2,...,A (120) 



where the 2N bosons <pj with integer and half-integer indices have the usual commutation 
relations [ipi(x), <fj(x')] = i7iSijSgn(x — x'). Clearly, Eq. ( |120|) is not a canonical transfor- 



mation because the number of degrees of freedom has been doubled. This formal difficulty 
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can be circumvented if we introduce an independent set of auxiliary bosons with half-integer 
indices Mi+1/2, i = 1,2, . . . , N with identical commutation relations, 

[u i+ i/ 2 {x), u j+ y 2 {x')) = in Kij sgn(x - x'), [ui{x), u j+ i/ 2 (x')) = 0. (121) 

Then, we can further write for the 131 state 

Uj+1/2 = ~<Pj + ¥3+1/2 ~ tPj+i, (122) 

and the transformation Uj,Uj + i/2 1— > fj becomes canonical. 

At the end of a calculation the auxiliary degrees of freedom need to be projected out. Such 
a projection has been explicitly carried out by us in Ref. jD]] for the simpler case of correlated 
quantum Hall bilayers. However, since the auxiliary degrees of freedom are independent 
of the physical ones, all quantum-mechanical averages involving physical quantities will 
automatically be correct, independent of the Hamiltonian chosen for the additional bosons 
Ui+i/2- In the following we select this Hamiltonian in the form (|20|) , i.e., identical to that for 
the physical bosonic fields sans the tunneling term. In this case, the transformations ( |1 20| ) 
and ([T22D give 

fL/2 ( 1 n _ 
Hi3i = dx I — v ij {'■ dx^i dxPj ■ + ■ d x (p i+ i/2 d x f j+1 / 2 : 



+E 



(27ray 



: e 



K<Pj-l/2 + ¥j - ¥>j+i ~ ^j+3/2) . + h c 



where the modified interaction matrix is given by the matrix product Vij = KnVij. Now we 
can use the standard fermionization prescription: 

^•+ a = -i=e^, (123) 
V zna 

where a = 0, 1/2. In particular, in the fermionic representation the original electron anni- 
hilation operator is written as a product of three fermion operators, oc ipj_i/2il ! j' l l ! j+i/2, 
while the complete fermionized Hamiltonian for the surface of the 131 system becomes 

rL/2 ( 1 



^131 = / dx(-J2 V n '■ *Pj+Jd x ip j+a : +- % ■ ^l+a^i+a ■ ■ ip} +a ip j+Q : 

+ E [A ^-1/2^^+1^+3/2 + h.c.]}. (124) 



A curious feature of this transformation is that the introduced auxiliary degrees of freedom 
can be thought of as located at the layers sandwiched between the physical layers — this is 
the reason for assigning half-integer indices to them. 

Similar fermion representations also exist for other correlated multilayer states. For the 
bosonic state 020, the interlayer correlations are absent, and the commutation relations can 
be diagonalized by the transformation 



u 



j = <^'T + <PjU ( 125 ) 
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where the auxiliary degrees of freedom labeled by the pseudospin index are located in the 
same layers. The fundamental boson operator can be written as tyj = ip^ipji, and the full 
Hamiltonian is 



^020 = / dx {- Vtt '■ Tpjjdxlpjv -+\ Yl ^iJ>ia ■ ■ Ipja^ja ■ 

+ E [ A ^Vll ^+u^+it + h-c] }, a =|, I . (126) 



Despite the similarity with its counterpart for the 131 state [Eq. (|124|) 1, this Hamiltonian 
involves subtly different correlations. 

In the case of 020 state, one can also avoid doubling the number of degrees of freedom 
by using the alternative fermionization introduced in Ref . . Specifically, the neighboring 
layers are combined into pairs, 

U2n+1 = ^2n+l + V 9 2n+2, U2n+2 = y?2n+l — ^2n+2, (127) 

and the tunneling operators are fermionized alternatingly as anomalous two- or four-fermion 
combinations, 

—L- e ^-^) = 4id x 4, e**-«»> = 4^s^4, • • • , (128) 



where ipj is given in Eq. ( |123| ). In cases where the number of layers is small, e.g., 
N = 2,3,4, alternative fermionizations exists for the 020 and 131 states, which in some 
cases allow exact solutions, see Ref. JTTJ] and Appendix [B]. 

The fermionization is also very simple for the pyramid states. As an example consider 
the 12321 state where the commutation relations are satisfied if we write 

Uj = ipj-i + (fj + ipj+i, (129) 

and the corresponding tunneling operator is bilinear, 

J_ e i(«,-i*n) = i/, J+1 . (130) 
zna J 

The edge theory with tunneling is quadratic and therefore exactly soluble, except for the 
complication of the soft modes mentioned in Section [II. 



VII. SUMMARY 

We have investigated the fractional quantum Hall effect in infinite layer systems. In the 
bulk we find quasiparticles with several unusual (irrational) features and at the edge we find 
surface phases which can be considered "chiral semi-metals" when the interlayer tunneling 
is not relevant. In addition we have presented fermionizations of the edge theory which may 
prove useful in finding exact solutions. 
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APPENDIX A: MODEL PROBLEM WITH MARGINAL TUNNELING 



Here we consider a model problem with marginal interlayer tunneling. The Hamiltonian 
is chosen in analogy with the Hamiltonian ([74]) , 



n 



^ dx [wpjid^j + X - [A e^ {x) (VH^'+i - id x ip} + h. 



(Al) 



with the marginal tunneling operator in parenthesis of a form reminiscent of the two-fermion- 
fermionized version of the marginal tunneling operator ( |128| ). The great advantage of the 
Hamiltonian ( |A~T| ) is its linearity, which allows us to calculate the conductivity for this 



model exactly. In order to do this, we perform a gauge transformation ipj 
Q + i — (j = jj. This gives, in obvious matrix notation, 



n 



L/2 
-L/2 



dx 



^Aid x ^ + -¥ (AV + VA) V 



<j( x )ip- with 



(A2) 



where ^ = (-^l ip2 ■ ■ -) T , the tri-diagonal matrix Aw = v5w + A(<5^/ + i + Si+iji) is coordinate- 
independent, while V = diag(Vi, V2, . . .) with Vj = d x Q. The corresponding transverse 
current operator has the form 



L/2 
-L/2 



dx^ ] 



iid x + - (tV + Vi 



(A3) 



where Uj = iX (S iJ+1 - 6 i+1)j ). 

The formal solution of this model (for a given realization of disorder) can be written in 
the limit of infinite system size (L — > 00) using the transfer matrix approach described in 



Appendix A of Ref. [12]. Using the Kubo formula for the vertical conductance, we obtain 



1 f L / 2 r~ 

I . ,_dx {n k - n k+UJ ) Tr V(x)S k (x, y)V(y)S k+UJ (y, x) 



2Nu 



L/2 



where n k = [eP vk + 1) 1 J the matrix in the vertex 

1 



V(x) = A- 1/2 
and the transfer matrix 



tid x + - (tV + Vt 



A" 1 / 2 , 



(A4) 



(A5) 
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S k (a,b) 



T x exp ^—i J dx 



kA' 1 - l -A- 1 ' 2 (tV + Vt) A' 1 ' 2 



(A6) 



To compute the disorder averages we assume that Vj(x) are independent delta-correlated 
Gaussian random variables with zero average, Vi{x) Vj(y) = A5ij5(x—y). Then, a calculation 
at a finite temperature T gives 



G z 



L 
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u + ^t 2 ) r^m^m 



(A7) 



where A(g) = v + 2Acosg, t(q) = 2Asing, and 



A(q) 



A 

T 



2A{q) + 1 



A 2 (q) 
Vv 2 - A 2 



As in the physical marginal cases, 020 and 131, at zero temperature and in the absence 
of disorder the conductivity vanishes. In the leading order in A Eq. ( |A7|) is identical to 
the result for 020 with delta-correlated disorder (0). However, in the clean limit we find 
G zz oc S(uj) for the model problem, which is not the correct behavior for the 020 and 131 
states. The difference arises because, in contrast to 020 and 131, in the model problem the 
current operator I z (|A3|) commutes with the full Hamiltonian (|A2|). 



APPENDIX B: EXACT SOLUTIONS FOR N = 4 

In this appendix we describe exact solutions for the spectrum of the edge theory for the 
special case of = 4 layers with periodic boundary conditions. We begin with a discussion 
of the 010 state with both tunneling and interactions and then consider the 020 state. We 
remark that similar solutions exist for the 010, 020, and 131 states in A^ = 3 layers. 

1. 010 State 



The Hamiltonian for the 010 state is 



K 



N=4 
010 



L/2 
-L/2 



dx 



:(d x Ui) 2 : +— :d x Uid x u l+1 : +— f e -<(«<-««+i) + h.c. 
4vr v ; 2tt + 2vra v 



^Bl) 



where u?, = u\. A straightforward fermionization ipi = e iUi ^ x ' / \ / 2na would yield a quartic 
Hamiltonian because of the density-density interaction term. We thus perform the orthog- 
onal transformation Ui(x) = Aij<pj(x) where 



A 



fl 1 1 1\ 

1 1 -1-1 
1-11-1 

\i -1-1 1 j 



;b2) 



The Hamiltonian ( |Bl"| ) then reads 
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^010 = / dx 



L/2 

A 



— Vi : (d x (pi) 2 : 

47T 



_|_ __ ( e -*(02-<fe) _|_ e -«(03-04) _|_ g«(</>2+03) _|_ g-*(03+</>4) _|_ n c 

27ra 



(B3) 



where t>i,3 = t> ± 2g and t> 2 = v 4 = t>. The original fields, Ui, are all radius R = 1 chiral 
bosons. We restrict ourselves to the sector where the total topological charge of the fields 
Ui is constant, which corresponds to the conservation of the total fermion number in the 
fermionic form of (|B1| ). Then the requirement that the topological charges of the original 
(u{) and transformed (4>i) bosons be integral is consistent with taking the fields (pi to all 
have unit radius as well. We can therefore fermionize according to rji{x) = e 1 ^^ /y/2na and 
then express the fermions in terms of their Fourier components f]i(x) = J2k etkx °ik/ VL. In 
terms of these mode operators the Hamiltonian (|B3j) is 



^010* = H [vikc\ k c ik + A (c\ k c 3k + c\ k c ik + c 2k c 3 -k + c\ k c\_ k + h.c.)] , (B4) 
k 

where the momentum k takes values in (27r/L)(Z + 1/2). This form of the edge theory is 
quadratic and hence exactly soluble via a Bogoliubov transformation. We find: 

<o 4 = E^fe, (B5) 
k 

where {b ik , b^ k ,} = Sij5 kk ' are four independent branches of fermionic oscillators with disper- 
sions 

ei(k) = {v + 2g)k, e s {k) = (v-g)k- } /{gk)* + W, 
e 2 (k)=vk, e 4 (fc) = (v-g)k+ ^(gk) 2 + 4\* 

We find that two of the branches develop curvature if and only if both g and A are nonzero. 



2. 020 State 



Now consider the 020 state whose bosonic Hamiltonian is given above in Eq. (p3[ ). The 
fields Ui all have unit radius, but their commutators are not conventionally normalized, 
see Eq. (pip . From our solution above for the 010 state we know that we can perform an 
orthogonal transformation ( |B2[ ) and still have four radius one bosons (pi. If we then define 
rescaled fields 9 Ax) = 6Ax)7y2 the 020 Hamiltonian becomes 



K 



N=4 
020 



L/2 



dx 



L/2 



ill 



:{d x 9i 



+ 



(2vra) ; 



e -iV2~(02-o 3 ) _|_ e -iV5(e 3 -6 4 ) _|_ e n/2(0 2 +0 3 ) _|_ e -iV2(e 3 +e A ) + ^ c 



(B7) 



where 9i are radius R = 1/ y/2 chiral bosons with conventionally normalized commutators, 
and the velocities Vi are the same as those given above in the solution of the 010 state. We 
can define a quartet of su(2)% KM currents as in Eq. Q54]). The Hamiltonian becomes 
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njN=A 
rl 020 



L/2 
-L/2 



dx 



1TV 



^ :[J t (x)] 2 : +4\(J*J* + r 3 jp j 



(B8) 



where we have again used the identity flSEp. Next, we define a new set of currents J°/{x) = 
R ab J b (x), which also obey independent su(2) 1 algebras, where R ab = S az 5 bx + 5 ay S by — 5 ax S bz . 
This rotation leaves the first term in the Hamiltonian ( |B8| ) unchanged and in the second 
term gives Jf \— > — J? . If we express the rotated KM currents, J?, in terms of four new 
radius R = 1 / -\/2 chiral bosons Oi via Eq. ([54]) , we find 



AT=4 
020 



i/2 1 ~ ~ 

rfx —Ma : d x 8id x 0i : , 

-L/2 47T J 



(B9) 



where 



M = 



fv/2 + g ^ 

u/2 A/tt 

A/tt u/2-p A/tt 

V A/tt u/2 / 



(BIO) 



We have succeeded in finding a quadratic representation of the Hamiltonian which can be 
readily diagonalized by performing an orthogonal transformation from the fields 6i to new 
fields dj. The final form of the Hamiltonian is 



n 



N=4 
020 



L/2 1 

dx—Vi :(d x di) : 

-L/2 47T 



(Bll) 



where the velocities are 



v 1 = v/2 + g v 3 = v/2-g/2 + . ] /(g/2y + 2(\/ny 
v/2 v 4 = v/2 - g/2 - Mg/2f + 2(A/tt) 1 



^2 



(B12) 



The exact spectrum of the 020 state contains four independent free chiral bosons whose 
velocities are renormalized by both interactions and tunneling. 



APPENDIX C: KLEIN FACTORS 

In this appendix we discuss Klein factors which are needed to produce the proper anti- 
commutation relations between different species of fermions. One place where Klein factors 
are needed is in the definition of the physical electron operators at the edges of each layer. 
We modify the definition of the electron operator at the edge of layer % ( p2"D to read 

Vi(x) oc e *^^ e i«i(*) j (CI) 

where A is an iV x iV matrix and Nj is the topological charge of the Bose field Uj(x) defined 
by 

Nj = -!- dxd x Uj(x). (C2) 
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Using the commutation relations of the chiral bosons (|21"D one can readily show 



and the combination appearing in the tunneling operator (34) can be written 

If we write A = nBK~ l , for some matrix B, then we have 

^(x)^(x') = Vjix'^iix^^e-™^- 3 ^, (C5) 
oc e i ["*( x )- u *+i(^)] e i7r ( B l +i,i+i- B l +i,i) e i7 i"( B i!- B i+i,!)(-^~ 1 )!fc Ar fe_ (qj) 



From Eq. ( |C5| ) we see that for the electron operators in different layers to anticommute 
we must have (B^ — Bji) e Z even if is odd and [B^ — Bji) e Z odd if if^- is even. From 
Eq. ( |C6| ) we see that if the elements of the matrix B are integers and (Bu — 
is an integer for all i and k, then the tunneling operator will be the same as in the absence 
of the Klein factors up to a sign which depends on the topological charge sector. 

As a concrete example, consider the multilayer 131 state. If the number of layers N is 
such that det(K) is odd (i.e., iVmod3 7^ 2) then the choice 

B fdet(*r) forj-z>2 
13 \ otherwise 

gives {^i(x), tyj(x')} = for all i,j and modifies the tunneling term by a factor of ±1. Note 
that in correlation functions, such as the current two-point function used to compute the 
conductivity ([77]), the tunneling operator is always accompanied by its hermitian conjugate 
and therefore factors such as e tn( - Bil ~ Bl+1 - l ^ K "> lkNk will cancel. 

In the above discussion we have constructed Klein factors out of the topological charge 
operators of the Bose fields in order to give the correct anticommutation relations between 
the physical electron operators without significantly modifying the tunneling Hamiltonian. 
When considering the fermionization of the multilayer edge theory, as discussed in Sec- 
tion 0, we have a different goal. In addition to ensuring that the electron operators an- 
ticommute we want the different species in the fermionized Hamiltonian to obey proper 
anticommutation relations. In this case an alternative approach to the one described above 
proves advantageous. This procedure is best described by considering a specific case, for 
example the 131 state. One can introduce 2N unitary operators F i+Q , where a = 0, 1/2 and 
i = 1, . . . , N, by formally enlarging the Hilbert space. These operators commute with the 
Bose fields Ui and obey 

{F i+a , F j+a ,} = {F i+a , F j+a ,} = {F i+a , F j+a ,} = 0, (C8) 

for % + a 7^ j + a'. The electron operators are modified according to 

%(x) oc F i _ 1 / 2 F i F i+ i/2e iu ^ x \ (C9) 



It is a straightforward exercise to demonstrate that Eqns. flC8|) and ( |C9|) imply that electron 
operators in different layers properly anticommute. The tunneling operator is now 
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oc ^- u ^F 3 _ V2 F 3 F} +1 F} +z/2 . (CIO) 

If we then follow the procedure described in Section [VT| with the modification that the new 
fermion species each absorb a factor of F: 

4 j+a = -t= Fj+a e^ x \ (Cll) 
V lira 

we arrive at the same fermionized Hamiltonian ( |124j ), but now with the different fermion 
species properly anticommuting. 



APPENDIX D: PROOF THAT Uea zz (u) oc 5(w) «-> [I Z ,H] = 

In this appendix we prove that the real part of the z-axis conductivity is proportional 
to a delta function in frequency if and only if the current operator commutes with the full 
Hamiltonian of the edge theory. The basic connection is of course the Kubo formula (jail): 

1 1 poo 

*'-M = toinzL***V'®- IM W'>- (D1) 

One of the directions is trivial: if [I Z ,H] = 0, then I z (t) is independent of time, and hence 
the commutator in Eq. ( |L)1[ ) vanishes identically implying that for all uo ^ 0, 3?e a zz (u) = 
and therefore ^tea zz oc 5(u). 

We now prove the converse, i.e., we assume 3tea zz oc S(u). Inverting the Fourier trans- 
form in Eq. ( |D1| ) we find 

/oo 
due- iuJt cu^ea zz (uj). (D2) 

Using our assumption in this equation gives 

<[/W(0)]> = J2e-^(n\[P(t),F(0)]\n) = 0, (D3) 

n 

where \n) and E n are eigenstates and eigenvalues of H. Note that the sum in Eq. (|D^) is a 
linear combination of terms (n\[I z (t), l z (0)]\n) with non-negative coefficients which depend 
on p. Since the sum must vanish for arbitrary (3 we conclude that 

(n\[I z (t),I z (0)]\n) = 0, Vn. (D4) 

Writing I z (t) = exp(iHt)I z (0) exp(-iHt) and inserting a summation over a complete set of 
states Eq. (|D4j ) is equivalent to 

^sm[(E n -E m )t]\{n\I z (0)\m)\ 2 = 0, Vn, Vt (D5) 

m 

The quantities |(n|/ 5: (0)|m)| 2 are clearly non-negative. Since the sum in Eq. ( p5| ) must 
vanish for all n and at all times t we can conclude 

|(n|/ z (0)|m)| 2 = 0, Vn^m. (D6) 

This equation says that the current has no off-diagonal matrix elements in the basis of 
energy eigenstates, i.e., the current operator V is diagonal in the basis of eigenstates of H. 
Since this means J z and H can be simultaneously diagonalized it implies that they commute: 
[I z (0),H] = 0, completing the proof. 
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